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ABSTRACT 

Multiple star systems are commonly assumed to form coevally; they thus provide the anchor for 
most calibrations of stellar evolutionary models. In this paper we study the binary population of 
the Taurus- Auriga association, using the component positions in an HR diagram in order to quantify 
the frequency and degree of coevality in young binary systems. After identifying and rejecting the 
systems that are known to be affected by systematic errors (due to further multiplicity or obscuration 
by circumstellar material), we find that the relative binary ages, |Alogr|, have an overall dispersion 
cr |Aiogr|~0.40 dex. Random pairs of Taurus members are coeval only to within CT|Aiogr|~0.58 dex, 
indicating that Taurus binaries are indeed more coeval than the association as a whole. However, 
the distribution of |Alogr| suggests two populations, with ^2/3 of the sample appearing coeval to 
within the errors (eriAiogrl^O.lG dex) and the other ~l/3 distributed in an extended tail reaching 
|Alogr|^0.4-0.9 dex. To explain the finding of a multi-peaked distribution, we suggest that the tail 
of the differential age distribution includes unrecognized hierarchical multiples, stars seen in scattered 
light, or stars with disk contamination; additional followup is required to rule out or correct for these 
explanations. The relative coevality of binary systems does not depend significantly on the system 
mass, mass ratio, or separation. Indeed, any pair of Taurus members wider than ~10' (~0.7 pc) 
shows the full age spread of the association. 

Subject headings: stars: binaries: general, stars: evolution, stars: fundamental parameters, stars: 
Hertzsprung-Russell diagram, stars: pre-main sequence, stars: open clusters and 
associations : individual ( Taurus- Auriga) 



1. INTRODUCTION 

Stellar evolutionary models are critical for interpret- 
ing astronomical observations, but they are not well- 
calibrated for pre-main sequence (PMS) stars. Such cal- 
ibration requires the measurement of some or all of the 
fundamental stellar properties: age, mass, radius, lumi- 
nosity, and effective temperature. Ages are notoriously 
difficult to estimate (Mamajek & Hillenbrand 2008; Hil- 
lenbrand 2009), though they can be inferred indirectly 
from a membership in a stellar population for which the 
mean age can be determined. Stellar masses and/or radii 
require orbital monitoring of eclipsing or visual binary 
systems. The known pre-main sequence eclipsing binary 
systems sparsely sample parameter space due to their 
extreme rarity (e.g. Irwin et al. 2007; Stassun et al. 
2007; Stempels et al. 2008). Most PMS visual binaries 
have only partial orbits because young stars are distant 
and any systems which can be spatially resolved neces- 
sarily have wide separations and corresponding long pe- 
riods (Steffen et al. 2001; Duchene et al. 2006), though 
a handful of short-period systems are bright enough for 
interferometric techniques to be feasible (Boden et al. 
2005; Schaefer et al. 2008). In contrast to ages, masses, 
and radii, the luminosities and temperatures of stars are 
straightforward to infer from single-epoch observations, 
so they offer the best near-term prospects for systematic 
calibration of stellar models. The procedure can also be 
inverted: given a star's luminosity and temperature, a 
theoretical model can be used to estimate its age and 
mass, plus its radius can be estimated directly from the 
Stefan-Boltzmann law. 

1 California Institute of Technology, Department of Astro- 
physics, MC 105-24, Pasadena, CA 91125 



The standard procedure for calibrating models with lu- 
minosities and temperatures is to place two or more nom- 
inally coeval stars on an HR diagram. These stars should 
trace an empirical isochrone sequence, and this sequence 
can be compared to theoretical isochrones in order to 
test their consistency with observations. HR diagram 
analysis has provided many crucial insights into mod- 
els of stellar interiors, atmospheres, and evolution (e.g. 
White et al. 1999; Luhman et al. 2003). There are many 
systematic astrophysical effects that can complicate this 
analysis, including unresolved multiplicity, obscuration 
from circumstellar material (i.e. an envelope or edge-on 
disk), and veiling from accretion (at blue optical wave- 
lengths) or circumstellar disk emission (at near-infrared 
wavelengths) . These effects can yield ages with errors of 
an order of magnitude or more, so the samples used in 
this analysis must be inspected closely to identify stars 
that might be affected. Additional physics, such as stel- 
lar activity or convection efficiency, could also play a role 
(Chabrier et al. 2007; Stassun et al. 2008). 

The inverse procedure plays an important role in the 
study of binary systems (e.g. Hartigan et al. 1994; White 
& Ghez 2001; Hartigan & Kenyon 2003). It is commonly 
assumed that multiple star formation proceeds almost 
simultaneously, such that all stars in a bound multiple 
system are coeval. However, this assumption can be 
tested only by using the evolutionary models that re- 
quire calibration, so any apparent disagreement between 
binary component ages could be due to non-coevality or 
errors in the models. Emerging evidence from several 
young eclipsing binary systems shows that their compo- 
nent properties are inconsistent with any single age pre- 
dicted by pre-main sequence evolutionary models (e.g. 
Stassun et al. 2007, 2008), but the frequency and degree 
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of implied noncoevality is still unclear. 

In this paper, we estimate the ages of a large sample 
of stringently vetted young binary systems in the Taurus 
association (r ~l-2 Myr; d ~145 pc) in order to test the 
system components' relative coevality and the validity of 
theoretical isochrones in matching empirical HR diagram 
sequences. In Section 2, we describe our sample of binary 
systems, and in Section 3, we describe the stellar models 
and analysis techniques used to infer stellar ages. In 
Section 4, we show an HR diagram with all our sample 
members and identify likely contaminants. In Section 5, 
we test the coevality of young binary systems by adopting 
model-predicted ages. Finally, in two appendices, we test 
the evolutionary models using a large sample of likely 
single stars and using a subset of high-order multiple 
systems. 

2. THE SAMPLE 

The accurate determination of stellar parameters re- 
quires spatially- and/or spectrally-resolved observations 
that are not polluted by light from companions; our sam- 
ple is comprised of all known Taurus binary systems that 
have spectral types and flux ratios for at least two compo- 
nents. At the distance of Taurus (~145 pc), most binary 
systems with separations of < 100-200 AU can not be 
resolved in seeing-limited observations with typical reso- 
lutions of Given this limit, most sample members 
have been observed with HST, ground-based adaptive 
optics, or echelle spectrographs (distinguishing the com- 
ponents' spectra and relative fluxes). High-order multi- 
plicity can give the appearance of non-coevality, so we 
immediately omitted any binary pairs which contained 
an additional component without its own measured spec- 
tral type and flux. As we describe in Section 4, we also 
cut additional targets from this sample if they appeared 
to be affected by other systematic uncertainties. We list 
the observed properties of our sample members in Table 
1, along with the references used to infer spectral types, 
extinctions, and fluxes. Given the ad hoc nature of our 
sample, it is not complete and may be subject to biases, 
but the goal of our initial sample selection was to be as 
inclusive as possible. 

The wide components that could be resolved in seeing- 
limited observations were drawn from the sample stud- 
ied in our previous wide multiplicity survey (Kraus & 
Hillenbrand 2007a, 2009). Most of these stars had al- 
ready been identified as Taurus members, so we drew 
their observed properties from previous work by Kenyon 
& Hartmann (1995), Duchene et al. (1999), White & 
Basri (2003), White & Hillenbrand (2004), and Luhman 
(2004, 2006). All of these authors reported spectral types 
and extinctions, and we adopted NIR magnitudes from 
the 2MASS Point Source Catalog (Skrutskie et al. 2006) 
or from our own PSF fitting photometry of 2MASS atlas 
images (Kraus & Hillenbrand 2007a). The 2MASS pho- 
tometry for RW Aur AB and most of the Duchene et al. 
sample was unreliable since the system separations fell 
near the 2MASS resolution limits, so we adopted total K 
fluxes from 2MASS and the K band flux ratios reported 
by White & Ghez (2001) or Corrcia et al. (2006). We 
adopted an outer separation limit of 30" for identifying 
binary pairs. The binary population probably dominates 
among all pairs extending out to ~2' (Kraus & Hillen- 
brand 2008), but the frequency of chance alignments be- 



tween unbound stars becomes significant at >30". 

The closer systems that could be resolved with high- 
resolution imaging were drawn from several recent spec- 
troscopic surveys. Most of the spectroscopic observations 
were obtained with HST/STIS by Hartigan & Kenyon 
(2003), but several systems were observed under good 
seeing by Duchene et al. (1999). Individual systems were 
also studied by White et al. (1999) with HST/FOS, or in 
our own survey of low-mass multiplicity with Keck laser 
guide star AO (LGSAO; Kraus et al., in prep). As be- 
fore, we inferred NIR magnitudes from the total system 
fluxes reported in 2MASS and spatially resolved flux ra- 
tios reported by Leinert et al. (1993), White & Ghez 
(2001), Correia et al. (2006), and our LGSAO survey. 

There are a small number of double- lined spectroscopic 
binaries that have been identified and studied in some 
detail in Taurus. We have added the well-known sys- 
tem UZ Tau Aab to our sample, adopting the spectral 
types and inferred H band flux ratio found by Prato et 
al. (2002). We also added the short-period spectroscopic 
binary V773 Tau Aab, which was studied with RV and 
intcrfcrometric monitoring by Boden et al. (2007) as 
part of an orbit monitoring program. The V773 Tau 
system also includes two faint companions at small sep- 
arations and a wide brown dwarf companion; we do not 
include the close companions because they do not have 
spectral type determinations and because they are too 
faint to influence the observed properties of V773 Tau 
Aab (Duchene et al. 2003), but we include the wide sub- 
stellar companion (Luhman et al. 2004). We also con- 
sidered whether to include V826 Tau (Massarotti et al. 
2005) and DQ Tau (Mathieu et al. 1997), but the only 
known flux ratio for V826 Tau is in the optical and there 
are no flux ratios reported for DQ Tau, so these systems 
could not be integrated with the rest of our sample. 

We also note that one binary system in our sample, XZ 
Tau, was recently suggested to be a possible hierarchical 
triple based on resolved radio observations of XZ Tau 
B at 1.3 cm and 7 mm. We have retained this system 
in our sample because the companion was not identified 
with HST (Krist et al. 2008) or in other surveys, but the 
suggested separation (^90 mas or 13 AU) was near the 
detection limits of the earlier observations, so we regard 
it as a credible possibility. As we will show later, XZ 
Tau B is overluminous compared to XZ Tau A or HL 
Tau, which supports its probable multiplicity. 

Finally, most of the binary components in our sam- 
ple have spatially resolved photometry in the K filter 
only. Inferring the component luminosities from such a 
red bandpass might introduce systematic errors in our es- 
timated luminosities due to near-infrared excesses from 
circumstellar disks. In order to address this prospect, 
we have searched the literature to determine which stars 
are likely to host a warm disk; we summarize our as- 
sessments and the corresponding references in Table 1. 
We based these assessments, in order of priority, on 3- 
10 /im photometry (from Spitzer/IRAC or ground-based 
AO imaging), optical spectroscopic accretion signatures, 
10-30 /im spectroscopy (from Spitzer/IRS), and finally 
on submm/mm photometry In each case where suffi- 
cient data is available, we have concluded that the star 
either has a disk ("Y"), does not have a disk ("N"), or 
might have a disk (but the observations aren't spatially 
resolved, so we can't determine which binary compo- 
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nent(s) have one; "Y?"). We will address the significance 
of flux excesses from warm dust in Section 5. 

3. ANALYSIS 

3.1. Inferred and Calculated Stellar Properties 

Any comparison of observations to stellar evolution- 
ary models requires the conversion of observed prop- 
erties (spectral types, filtered magnitudes, and extinc- 
tions) into fundamental physical parameters (effective 
temperatures and bolomctric luminosities). This pro- 
cess is accomplished by invoking a temperature scale, 
a set of bolometric corrections, a reddening law, and 
an estimated distance. Temperature scales directly re- 
late spectral types to temperatures, allowing the esti- 
mation of temperatures based on observed spectral fea- 
tures. Temperatures are also used to define intrinsic col- 
ors, from which an observed color can be used to in- 
fer the reddening and extinction. Bolometric corrections 
are temperature-dependent ratios of the flux in a filtered 
band to the full bolometric flux and are calibrated for 
nearby field stars that have been studied across the full 
range of wavelengths with significant contribution to the 
luminosity. Once the bolometric flux is known, the dis- 
tance for a star then directly yields the bolometric lumi- 
nosity. 

Observed spectral types can be converted into effec- 
tive temperatures using an adopted temperature scale, 
but this process is somewhat uncertain for young stars 
since the relation between temperature and spectral type 
is sensitive to surface gravity. For M-type stars, a giant 
with a given spectral type can be ~300 K warmer than 
a dwarf of the same type (e.g. Leggett ct al. 1996 versus 
Perrin et al. 1998; Richichi et al. 1998; van Belle et al. 
1999). Young stars have intermediate surface gravity, so 
it has been suggested (e.g. Luhman 1997; Luhman et al. 
2003) that an intermediate temperature scale might be 
appropriate. For example, Luhman et al. (2003) found 
that M2-M5 stars in IC 348 have bluer R-I colors than 
their field counterparts, matching the trend seen for gi- 
ants. Similarly, we have shown in a previous publication 
(Kraus et al. 2006) that the V - 1 colors of M5-M9 stars 
in Taurus and Upper Sco are significantly bluer than 
those of field dwarfs, though not as blue as for giants. 
Conversely, more recent surveys of moderately older stars 
in e Cha and v Cha by Lyo et al. (2004, 2008) show that 
the broadband colors and temperature-sensitive narrow- 
band spectral indices are generally similar to field dwarfs 
to spectral types as late as M5.5, with only a mild dis- 
crepancy in the B — V color that could be attributed to 
chromospheric activity. 

We have adopted the temperature scale suggested by 
Luhman et al. (2003) for use in low-gravity young stars. 
For spectral types <M0, the temperature-SpT relation 
does not appear to be gravity-sensitive, so Luhman et al. 
use the temperature scale of Schmidt-Kaler (1982). For 
late-type (>M0) young stars, Luhman et al. define an 
intermediate temperature scale that makes the average 
cluster sequences of Taurus and IC348 internally coeval 
with respect to the NcxtGcn models (Baraffc ct al. 1998). 
Unfortunately, this calibration introduces an element of 
circularity into our analysis; the temperature scale is cho- 
sen by definition to make our chosen models reproduce 
observations. It is unclear whether the temperature scale 
of young stars is truly different or the adopted temper- 



ature scale is correcting for a discrepancy in the model, 
so our results should be weighed accordingly. The dif- 
ference in inferred temperatures with respect to field M 
dwarf values (Leggett et al. 1996) is -50 K at M3, -100 
K at M5, and ~200 K at M8, equivalent to a systematic 
uncertainty of ±0.25, ±0.5, and ±0.75 subclass, respec- 
tively. 

We have adopted the bolomctric corrections that we 
previously described for use with field stars (Kraus & 
Hillenbrand 2007b) . For spectral types <K7, we used the 
corrections suggested by Masana et al. (2006), while for 
M dwarfs, we used the corrections of Leggett et al. (1992) 
and Leggett et al. (1996). Values calibrated for field 
dwarfs might be systematically inaccurate if the bolo- 
metric correction is also gravity-sensitive, so it would be 
helpful to verify whether these bolomctric corrections are 
valid for young stars. However, we are not aware of any 
such tests having been attempted. In all cases, we have 
adopted the observed NIR flux (and its bolometric cor- 
rection) which is closest to the J filter. As we describe 
below, luminosity excesses are typically least significant 
in the I and J filters, but most of the binary systems in 
our sample have resolved flux measurements only in the 
H or K filters. 

We have implemented our extinction corrections using 
the interstellar reddening law of Schlegel et al. (1998), 
which stipulates that one magnitude of visual extinction 
corresponds to Aj = 0.28, A H = 0.18, and A K = 0.11; 
these values are consistent with the interstellar redden- 
ing law in the 2MASS filters suggested by Indebetouw et 
al. (2005). Reddening laws might vary in regions with 
extremely high density (Weingartner & Draine 2001; 
Roman-Zuniga et al. 2007), but most of our sample 
members are only moderately reddened. Analysis of 
2MASS source counts and colors toward the Ophiuchus, 
Lupus, and Pipe Nebulae suggest that the standard in- 
terstellar reddening law is appropriate to extinctions of 
Ay^ 20 or more (Lombardi et al. 2008). 

We adopted a characteristic distance for all Taurus 
members of 145±15 pc. Recent high-precision parallax 
measurements with the VLBA (Lestrade et al. 1999; 
Loinard et al. 2007; Torres et al. 2007; Loinard et al. 
2008) suggest that there might be a distance gradient 
of 165-125 pc in the east-west direction, though the dis- 
crepant distances of neighboring V773 Tau and Hubble 4 
(148±5 pc versus 132.5±0.6 pc) suggest an overall scat- 
ter at any location of ~10-15 pc (10%). The luminosity 
uncertainty if we adopt the characteristic distance is only 
~0. 1-0.2 mag, which is similar to the uncertainty from 
dereddening and intrinsic variability, so attempting to 
extrapolating more precise distances from this sugges- 
tion of 3D structure is not likely to improve our results. 

We should note that there is room for significant un- 
certainty in our inferred luminosities due to the intrinsic 
variability of young stars. Variability in Class III stars 
should be caused by cool spots, so its characteristic am- 
plitude at near-infrared wavelengths should be no more 
than -0.1 mag (e.g. Carpenter et al. 2001, 2002). How- 
ever, the same survey showed that Class I-II stars occa- 
sionally vary by as much as 1-2 mag in the near infrared, 
and extreme classes of stars (such as FUor and EXor 
stars) can vary by even more. This suggests that an un- 
usually young apparent age for one binary component 
could be the result of variability in that component. 
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A star that is surrounded by circumstcllar material 
could also appear systematically undcrluminous by sev- 
eral magnitudes if it possesses an edge-on disk or a mas- 
sive circumstellar envelope that renders it visisble only 
in scattered light. It would then appear much older than 
its unobscured companion. One such case is HK Tau B; 
both components of this system have a spectral type of 
Ml (White & Hillenbrand 2004), but the B component 
is 3.3 magnitudes fainter in K. 

Finally, the presence of veiling due to accretion (at 
blue wavelengths) or circumstellar dust emission (at red 
wavelengths) could bias the inferred luminosity of a bi- 
nary component. Observations in the I or J filter are 
least sensitive to this veiling, though even those measure- 
ments can be significantly impacted (Folha & Emerson 
1999; White & Hillenbrand 2004). In our study, we have 
used photometry from J or from the nearest redder fil- 
ter. We would prefer to use J for all systems, but in 
most cases, the absence of suitable data forces us to use 
relative photometry in H or K . Data from bluer filters is 
also available for some sample members, but these data 
are not as homogeneous and we prefer to minimize the 
total wavelength range over which we make our luminos- 
ity estimates. 

We list the stellar properties that we inferred with 
these methods in Table 1. Our temperature uncertain- 
ties have been determined from the uncertainty in each 
star's spectral type. Each of the bolometric magnitude 
uncertainties listed above should contribute ~0.1 mag, 
so we have adopted a total statistical uncertainty of 0.3 
mag. However, many of the uncertainties (such as for 
distance and extinction) should be correlated between 
binary components, so the uncertainty in their relative 
ages should be lower than our formal estimates. As we 
discuss further in Section 5, there is compelling evidence 
that the scatter in relative ages of binary components is 
indeed significantly smaller, suggesting that these error 
estimates are conservative. 

3.2. Inferred Physical Stellar Parameters 

Several sets of pre-main sequence evolutionary models 
have been developed in recent decades, but all of these 
models still face significant challenges in confronting ob- 
servational constraints. Hillenbrand & White (2004) 
found that all models have difficulty matching the dy- 
namical masses of young stars with M <1.2 Mq, a range 
which encompasses almost all of our sample. However, 
the Lyon models (Baraffc et al. 1998; Chabrier et al. 
2000) seem to work best for low-mass stars, especially 
when using a mixing length of a= 1.0 and the revised 
temperature scale of Luhman et al. (2003). All models 
were found to reproduce the observed masses for stars 
with M >1.2 Mq, though the Lyon models only extend 
to 1.4 Mq. 

In light of these results, we have adopted a hybrid 
combination of the Lyon models for low-mass stars and 
the models of D'Antona and Mazzitelli (1997; DM97) 
for higher-mass stars. For masses <0.5 Mq, we use the 
mass-luminosity-temperature relations of the Lyon mod- 
els with a mixing length of a= 1.0. For masses >1.0 Mq, 
we use the corresponding relations of DM97. Finally, 
in the intermediate regime of 0.6-0.9 Mq, we adopt a 
weighted average of the luminosity and temperature pre- 
dicted by each model in order to produce a smooth tran- 



sition between the two sets of models. This choice does 
not provide any insight into the missing physics that are 
still required to bring the models into agreement, but it 
represents an acceptable compromise for estimating rel- 
ative ages of a sample of young stars. The two sets of 
models converge at older ages, so our solution is only 
important for very young stars (<10 Myr). 

These models report stellar luminosities and tempera- 
tures at quantized values of age and mass, so for each of 
our sample members, we have linearly interpolated be- 
tween the four values of T e f / and Mboi around it in the 
HR diagram. The Lyon models also face a significant 
challenge with respect to very young stars since they are 
not defined for ages of <1 Myr, so for each star that falls 
above this isochrone and has a mass within the affected 
range (<1 Mq), we have linearly extrapolated its age 
from the four points below it in the HR diagram. These 
extrapolated ages should be regarded as much more un- 
certain than older ages, but the degree of error should 
be similar for stars with similar HR diagram positions, 
so only systems with disparate masses will be subject to 
the full systematic uncertainty. 

In Table 2, we list the inferred mass and age for each 
star from our hybrid system and from the default Lyon 
and DM97 models. Several sample members illustrate 
the extreme difference in mass and age estimates for the 
two sets of models. For example, the inferred parame- 
ters of HN Tau A arc M =1.35 Mq and log(r) =6.85 
according to the Lyon models and M =0.65 Mq and 
log(r) =6.05 according to the DM97 models; our hybrid 
isochrones yield M =0.85 Mq and log(r) =6.27. The un- 
certainty in log(r) due to observational errors is ±0.25 
dex for an average star in our sample, but as we describe 
below, this uncertainty could be an overestimate since 
binary systems appear to be more coeval. 

4. THE HR DIAGRAM 

The overall population sequence in an HR diagram pro- 
vides a valuable test of our choice of evolutionary models. 
If the association is nominally coeval, then it should trace 
a single recognizable sequence that is parallel to theoreti- 
cal isochrones. Individual sample members that strongly 
deviate from the association sequence should also be ex- 
amined for a systematic source of error such as erroneous 
membership, circumstellar material that blocks and scat- 
ters the stellar flux, or misclassification of spectral types. 

In Figure 1, we show the HR diagram for all of our 
binary sample members. The cluster sequence seems to 
trace the 1-2 Myr isochrone, albeit with significant scat- 
ter due to the many sources of uncertainty and the un- 
known spread of stellar ages. We adopted the two sets 
of theoretical isochrones (Lyon for low-mass stars and 
DM97 for higher-mass stars) specifically because they 
fit the Taurus single-star sequence across the full mass 
range of our sample (Appendix A), so this agreement 
with the canonical age of Taurus is not surprising. How- 
ever, many individual members fall unusually low; as we 
show in Figure A2 for single stars, the 5 Myr isochrone 
represents a +2a deviation from the median age of Tau- 
rus, so we adopt this approximate lower limit for the 
observed scatter of the main body of members. Most of 
the binary components below this isochrone are known 
to be anomalous. 

The warmest anomalous members (HBC 352-357) have 
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been classified as Taurus members for several decades 
(e.g. Walter et al. 1988; Herbig & Bell 1988; Kenyon 
& Hartmann 1995) based on their X-ray emission and 
(for HBC 352-355) their radial velocities. However, these 
stars are located at the far western edge of Taurus, well 
away from the central cloud cores, and their undcrlu- 
minosity has been recognized since their discovery. Few 
membership surveys have extended so far in this direc- 
tion from the clouds, so it is unknown whether these stars 
are surrounded by a more extensive coeval population. 
Given their proximity in projection to the Perseus com- 
plex, it seems plausible that they are associated with that 
more distant, but similarly young population. Since both 
members of each binary pair seem to be equally anoma- 
lous, this seems to be a reasonable explanation for their 
low positions on the HR diagram; we therefore choose to 
remove them from our sample for all subsequent analysis. 

Three of the other binary companions arc known 
sources seen only in scattered light due to the presence of 
an edge-on circumstellar disk. Stapelfeldt et al. (1998) 
used HST and AO observations to show that the optical 
and NIR flux from HK Tau B comes from extended nebu- 
losity, with no recognizable flux coming directly from the 
central star. Krist et al. (1995) and Krist et al. (1998) 
found similar results from HST imaging of HL Tau and 
Haro 6-5B, respectively. HK Tau and Haro 6-5B appear 
significantly undcrluminous in our HR diagram since this 
reflected light only represents a small fraction of each 
star's total emitted flux. Surprisingly, HL Tau does not 
appear under luminous, which suggests that its luminos- 
ity might be dominated by scattered light in the opti- 
cal and direct flux from the central star in the J band. 
However, we have chosen conservatively to omit all three 
companions from our analysis of relative binary ages. 

The binary component Haro 6-28 A sits just below the 

5 Myr isochrone, so further analysis of its scattered-light 
properties might be worthwhile in the future. However, 
its companion Haro 6-28 B sits just above the 5 Myr 
isochrone, so the inferred ages of the two components 
are mutually consistent. Barring a systematic uncer- 
tainty for the binary system, this consistency suggests 
that Haro 6-28 AB is genuinely one of the oldest sys- 
tems in Taurus. Haro 6-28 was originally identified as an 
Ha emission line star (Haro et al. 1953; Cohen & Kuhi 
1979) and at least one component has a 1.3mm excess 
indicative of a disk, so its youth seems secure. 2 

Finally, two companions (V710 Tau C and 
2M04554801) sit below the 5 Myr isochrone with- 
out any obvious explanation. Their optical spectra 
appear to be accurately classified (Luhman 2004; Kraus 

6 Hillenbrand 2009), so a large error in temperature 
seems unlikely. They are associated with stars that 
appear youthful, so membership in a different population 
also does not explain their anomalously old apparent 
ages. V710 Tau C has not been observed at high spatial 
resolution, but it shows a very significant K band excess 
in 2MASS (J — K <~2.2), so it might possess an edge-on 
circumstellar disk. 2M04554801 appears to be a point 
source in K band imaging with Keck LGSAO (Kraus 

2 As was demonstrated by Carpenter et al. (2006) for Upper 
Sco, massive circumstellar disks arc relatively rare (/ ~5%) for 
~0.3-0.5 Mq stars by the age of ~5 Myr, though not as rare as for 
higher-mass stars. 



& Hillenbrand, in preparation), so if the star is seen in 
scattered light in the near-infrared, then the scattering 
region must be smaller than for other sources with 
edge-on disks. The original discovery spectrum shows 
obvious signatures of youth (K. Luhman, priv. comm.), 
and the optical/NIR SED does not show any of the 
characteristic signs of an edge-on disk, so the explana- 
tion for its undcrluminosity is currently unknown. We 
can not justify removing either star from our sample, 
but both systems should be regarded with appropriate 
skepticism. 

In the following analysis, we will omit all of the HBC 
sources that have questionable membership. We will re- 
tain Haro 6-5B and HL Tau for a more in-depth study of 
high-order multiple systems, but will omit all three con- 
firmed scattered light systems for testing coevality. We 
will retain Haro 6-28 A, V710 Tau C, and 2M04554801 
since there is no conclusive evidence to suggest that they 
suffer a systematic bias. 

5. THE COEVALITY OF YOUNG BINARY SYSTEMS 

5.1. The Relative Ages of Binary Systems 

The first step of our analysis is to determine whether 
binary systems appear more coeval than the association 
as a whole. The timescale for star formation across an 
entire region could be as long as 10 Myr (Mouschovias 
1976; Shu 1977), while the formation of a single star 
system should proceed on the dynamical timescale of 
-0.1-0.2 Myr (Shu et al. 1987). An upper limit on 
the non-coevality of binary systems will provide a direct 
constraint on the formation timescale for binary systems 
with respect to the entire association. 

In Figure 2 (top), we show the same HR diagram as 
Figure 1, minus the objects we eliminated in the discus- 
sion above, where each of the binary pairs is connected 
by a line. The overall trend for binary pairs is to de- 
fine lines that roughly parallel the theoretical isochrones, 
which is expected since the models have been chosen by 
comparison to the Taurus single-star sequence. However, 
as we determine quantitatively below, some pairs fall on 
significantly different isochrones. This type of plot pro- 
vides a summary of the underlying data, but it is hard to 
draw any firm conclusions regarding overall coevality or 
possible dependence of coevality on binary parameters. 
More detailed statistical analysis must be pursued using 
the inferred stellar ages, as has been concluded by past 
studies of binary ages (e.g Hartigan et al. 1994; White 
& Chez 2001; Hartigan & Kenyon 2003). 

In Figure 3 (top), we show a histogram of the abso- 
lute difference in the inferred logarithmic age, |A logr| = 
| log Tp r i m — log T sec , for each of our binary pairs. The 
RMS scatter in |Alogr| among our sample is 0.40 dex 
and should encompass all observational uncertainties as 
well as any intrinsic age spread for binary pairs. Most of 
the model-predicted masses fall in the range 0.3 < M < 
0.9 M , with some masses extending as high as 1.5 Mq 
and as low as 0.03 Mq. We have compared the scatter 
in binary component ages to that of the overall Taurus 
population by using a bootstrap Monte Carlo routine to 
simulate 10,000 populations where we pair each primary 
with another randomly-selected secondary. In Figure 3 
(bottom), we show the distribution of all RMS scatter 
measurements for these simulated pairs; only 14 realiza- 
tions of our simulation (~0.15%) have RMS scatter of 
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Fig. 1. — HR diagram for all components of all sample binaries. The binary components that we rejected due to known systematic errors 
(Section 4) are shown with open circles, while the rest of our sample is shown with filled circles. The dashed lines denote isochrones at 1 
Myr (red), 5 Myr (green), 10 Myr (blue), and 50 Myr (black). Most Taurus members fall along the 1 Myr isochrone, but 10 fall below the 
5 Myr isochrone. Three binary pairs that fall mutually below this limit might be associated with the more distant Perseus star-forming 
complex, while individual components that are associated with apparently young binary companions could be seen in scattered light or 
have incorrect spectral types. We use solid lines to connect each binary pair with one or more rejected components. 

0.40 dex or less, indicating that binaries are more co- 
eval than the overall Taurus population at a significance 
of ~3<r. The typical age differences for binary systems 
and for random pairs are similar to those measured by 
White & Ghez (2001); they used a similar sample, but 
estimated effective temperatures using dercddcncd V — I 
photometry instead of spectral types. 

As we describe in Section 3, the observational uncer- 
tainties for our sample allow us to estimate the corre- 
sponding uncertainty in each binary component's age, 
and thus in the degree of coevality. These estimated un- 
certainties vary significantly across our sample, but the 
median, mean, and quadratic mean of these uncertain- 
ties (0.33, 0.42, and 0.48 dex in measured |Alogr|, re- 
spectively) are all similar to the standard deviation for 
our sample. This strongly suggests that much of the to- 
tal error budget is dominated by observational errors, 
and therefore that model-related errors and the true dis- 
persion in relative ages for binary components are both 
<<0.40 dex. 

However, we must also consider whether a single dis- 
tribution is adequate for describing all binary systems, 
as there are several effects that could bias one binary 
component's age by a significant amount (including over- 
luminosity due to unresolved additional multiplicity or 



underluminosity because an object is seen in scattered 
light). We could expect a narrow distribution centered 
close to zero and broadened by the observational errors 
(corresponding to unaffected binary systems) plus a sec- 
ondary peak away from zero (for systems affected by un- 
usual phenomena such as those mentioned above). The 
observed distribution seems to match this expectation, 
with most systems concentrated at |Alogr|<0.3 dex and 
an extended tail at Alogr|>0.5 dex. If we omit the ex- 
tended tail and compute the standard deviation of the 
binary age dispersion for only systems with |Alogr|<0.4 
dex, we find a dispersion of C|Ai O gT|^0.16 dex, corre- 
sponding to a typical factor of 1.5 in relative age. This 
dispersion is actually lower than the typical uncertainties 
we estimated above for |Alogr|, which suggests that we 
might have been too conservative in estimating observa- 
tional uncertainties. 

The extended tail includes one sample member that we 
suggested to be a candidate edge-on disk (V710 Tau C) 
based on its underluminosity and extremely red J — K 
color. The tail also includes two possible hierarchical 
triple systems, DK Tau and XZ Tau, which were sug- 
gested to be possible hierarchical triples by Jensen et al. 
(2004) and Carrasco-Gonzalez et al. (2009). Polarization 
measurements by Jensen et al. (2004) indicate that the 
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component disks in DK Tau AB are misaligned, unlike 
most other double-disk systems in Taurus. Radio obser- 
vations by Carrasco-Gonzalez et al. (2009) show that XZ 
Tau B has a double-peaked distribution that could result 
from a ~13 AU binary companion. 

In light of this possible bimodality, it is worthwhile to 
return to Figure 2 and plot only the "coeval" sample 
(|Alogr|<0.4 dex; center) and the "non-coeval" sam- 
ple (|Alogr|>0.4 dex; right). This division reveals a 
startling trend; among the "non-coeval" subsample, 11 
of the 12 systems have a significantly younger primary 
star. As we show in Appendix A, the model-derived ages 
in this mass range do not show a mass-dependent trend, 
so the tendency for some binary primaries to appear 
younger must be either a genuine result of the formation 
process or a result of binary-specific systematic uncer- 
tainties. Given the clear discrepancy with respect to the 
apparently coeval majority of our sample, we strongly 
suspect that systematic errors are to blame. 

The fraction of apparently coeval binary systems in 
our sample (24/36) is identical to the fraction identified 
by the survey of Hartigan et al. (1994). In a sample 
of binary systems in Taurus and Orion (which included 
many Taurus binaries that we have rejected as hierarchi- 
cal multiples), they found that 17/26 had Alog(L)<0.24 
dex, corresponding roughly to |Alogr|<0.4 dex. On its 
face, this result suggests that all of the improvements 
in evolutionary models, spectral type assessments, and 
multiplicity surveys in the past 15 years have only served 
to cut the standard deviation in |Alogr| for coeval sys- 
tems from 0.23 dex to 0.16 dex. However, there is one 
significant difference. All of our non-coeval systems pos- 
sess apparently younger primaries, while all of the cor- 
responding systems in the Hartigan et al. sample pos- 
sess apparently older primaries. There is little overlap 
between our samples since many wide Taurus binaries 
have since been discovered to be hierarchical multiples, 
so one possible explanation is that our stringent multi- 
plicity vetting simply allows another systematic error to 
dominate. There are three likely culprits for the system- 
atic error that makes some systems appear non-coeval. 

First, the non-coevality could result from unrecognized 
high-order multiplicity. The binary fraction is higher 
among solar-type stars than lower-mass stars (Duqucn- 
noy & Mayor 1991; Fischer & Marcy 1992). If this trend 
also applies to the fragmentation of binary components 
into high-order multiples, then we might expect more bi- 
nary primaries to be unresolved pairs (which would then 
appear to be a single overluminous star). An ovcrlu- 
minosity by 0.75 mag (denoting an equal-mass binary 
pair) should correspond to an apparent age discrepancy 
of ~0.5 dex, which matches the observed trend. How- 
ever, our sample includes many objects with similar tem- 
peratures (and thus presumably similar masses), and as 
we will show in a future paper (Kraus & Hillenbrand, 
in preparation) , the frequency of binaries at separations 
<50-100 AU in Taurus is nearly constant for all masses 
>0.3 Mq. We might expect a significant excess of ap- 
parently younger primary stars if the secondary masses 
fell significantly below this limit, but few of our targets 
do. 

The other two possible explanations are tied to the 
properties of circumstellar disks. Surveys of protoplan- 
etary disks in binary systems have suggested that disks 



might be more likely to form or persist around the pri- 
mary than the secondary (e.g. Monin et al. 2007), 
even though lower-mass stars generally retain their disks 
longer (e.g. Carpenter et al. 2006). We determined most 
stellar luminosities from the K band flux, so a disk ex- 
cess could have led to significant overestimation of those 
luminosities. In addition, the disks of high-mass stars 
boast more substantial NIR excesses than those of low- 
mass stars and brown dwarfs (e.g. Meyer et al. 1997 
versus Liu et al. 2003), so the magnitude of the luminos- 
ity overestimate should also be larger for primaries than 
secondaries. 

All of these explanations should be investigated and 
ruled out before an astrophysical explanation is consid- 
ered. In particular, systematic effects from disks should 
be mitigated by estimating stellar luminosities using the 
least contaminated filter (J) and by modeling the circum- 
stellar dust emission using JHKL photometry so that 
any remaining excess can be subtracted. However, our 
preliminary disk census suggests that NIR excess con- 
tamination might play only a modest role in biasing rel- 
ative binary ages (unless the disk directly obscures the 
central star, making it appear in scattered light; Sec- 
tion 4). Of the 24 pairs which appear coeval, 20 have 
at least one disk (where >4 are mixed pairs and >11 
are double-disk systems) . This presents little contrast to 
the 12 apparently non-coeval pairs, of which nine have at 
least one disk (with 3 mixed pairs and 6 double-disk sys- 
tems). The similar and nontrivial fractions of mixed pairs 
are difficult to explain if disk biases dominate, though 
double-disk systems could appear coeval if both binary 
components' luminosities are biased upward by the same 
amount. 

5.2. The Role of System Parameters in Binary 
Coevality 

The detailed physics of multiple star formation are still 
poorly understood, so any apparent trends in the coeval- 
ity of binary systems could yield valuable new constraints 
on theoretical models. A third of the systems in our 
sample appear non-coeval, so any such trend could be 
identified among the 12 non-coeval binary pairs. The 
three binary properties that we can test against system 
coevality are the component mass ratio, the total system 
mass, and the system projected separation. The degree 
of coevality as a function of separation across the entire 
association could also constrain the large-scale star for- 
mation processes, so we also analyze the coevality as a 
function of separation between all pairs of stars in our 
sample. 

In Figure 4, we plot the difference in system age 
IAIoktI function of binary mass ratio. If binaries 
truly formed non-cocvally, then we might expect the sys- 
tems with the most extreme or most similar mass ratios 
to show the largest discrepancy in ages. However, the 
dispersion in |Alogr| for the five systems with q < 0.3 
is 0.41 dex, similar to the overall dispersion for our full 
sample (0.40 dex). If we limit this analysis to only the ap- 
parently coeval population (|AlogT|<0.40 dex), the dis- 
persions are 0.17 dex and 0.16 dex, respectively. This 
indicates that there is no strong trend for reduced co- 
evality in these extreme systems. Hartigan et al. (1994) 
also found no such trend in their sample. 

In Figure 5, we plot |Alogr| as a function of total 
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Fig. 2. — HR diagram for the binary pairs in our sample, where 
each pair is connected by a line. The top panel shows all systems, 
whereas following the text in Section 5.1, the other panels show 
only systems with |Alogr| < 0.4 dex (middle) and |AlogT| > 
0.4 dex (bottom). The binary systems in our sample trace the 
approximate contours of stellar evolutionary models, suggesting 
that the overall trend is correct, but our detailed results are more 
easily described in terms of the inferred stellar ages (Figures 3-6). 

system mass. Most of our sample spans only a limited 
mass range (0.7-1.5 M©), but we see no evidence of a 
mass-dependent trend. Our sample includes only four 
systems with a total mass of <0.5 M Q , but we also see 
no significant trend for a higher scatter in ages. The 
dispersion (0.36 dex) is similar to that of the full sam- 
ple, though almostly entirely dominated by one system 
(2M04554757+2M04554801). Our sample includes only 
two high-mass pairs, RW Aur AB and V773 Tau Aab, 
for which we measure age discrepancies of 0.55 dex and 
0.01 dex, respectively. 

In Figure 6, we plot |Alogr| as a function of system 
separation. If the separation of a binary system scales 
with the protostellar core size when fragmentation oc- 
curred, then wider systems should typically fragment at 
an earlier stage than closer systems. This suggests that 



wide pairs might fragment earlier and evolve more in- 
dependently, possibly yielding binary components with 
a larger dispersion in apparent ages. However, as for 
the previous figures, this comparison does not indicate 
any significant role of separation in establishing the bi- 
nary component ages. The inner and outer halves of the 
sample (divided at 800 AU) have dispersions of 0.37 dex 
versus 0.42 dex (for the full set) and 0.14 versus 0.18 dex 
(for the coeval subset). We conclude that binary sys- 
tems of all separations are similar coeval to within our 
observed limits. 

5.3. The Intra- Association Coevality of Young Stars 

The lack of a separation-dependent trend in differential 
age begs an important question. If binary pairs are sim- 
ilarly coeval with their associated components, but sig- 
nificantly more coeval than the association as a whole, 
then what is the form of the transition between these 
regimes? Are adjacent (but unassociated) Taurus mem- 
bers more coeval than distant members, or is the age 
spread similar across all spatial scales greater than the 
binary separation regime? The distribution of Taurus 
members has been suggested to represent a small num- 
ber of subclusters (e.g. Gomez et al. 1993) with radii 
of ~l-2° (~5 pc), though it is unclear whether those 
groups are distinct from the large-scale filamentary (and 
possibly fractal) structure (Kraus & Hillenbrand 2008). 
If these apparent groupings are closely associated, then 
stars separated by <5 pc might be more coeval than the 
wider association. 

In Figure 7, we address this question by plotting 
|Alogr| as a function of separation for all possible pairs 
of the primary and secondary stars in our binary sam- 
ple. We can not draw any conclusions for separations 
of 30-1000" (5000 AU to 0.7 pc) due to small number 
statistics, but the dispersion at larger separations (as in- 
dicated by the standard deviation in bins 0.5 dex wide) is 
consistently ^0.6 dex across the entire separation range. 
This result suggests that the coevality we see for binary 
systems (Figure 6) is limited to scales of <0.7 pc. We 
can not determine if this coevality is limited exclusively 
to binary systems, though; as we showed in our analy- 
sis of the spatial distributions of stars (using two-point 
correlation functions; Kraus & Hillenbrand 2008), the bi- 
nary regime only encompasses separations of <2'. The 
observed scatter is unlikely to result from any distance 
dispersion of Taurus members (<15 pc for members in 
similar parts of the cloud; Section 3.1) since it would 
yield a scatter of <0.2 mag in Mf, i or ^0.15 dex in logr 
(0.20-0.25 dex in |Alogr| for a pair of stars). 

5.4. Implications for (Multiple) Star Formation 

Our results for binary pair age differences are con- 
sistent with theoretical predictions for the timescale of 
local star formation (i.e. within one core) and global 
star formation (spanning the entire molecular cloud). 
If most binary systems in Taurus are coeval to within 
<0.16 dex (including observational uncertainties), then 
given its median age (1.8±0.2 Myr; Appendix A), the 
formation times for binary components must differ by 
<0.7 Myr. The expected timescale for an individ- 
ual protostar to collapse after achieving supercritical- 
ity is the dynamical timescale (^0.1-0.2 Myr; Shu et 
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Fig. 3. — Top: Distribution of differences in logarithmic age, |Alogr|, for all 36 pairs of stars in our sample. The RMS scatter in Alog-r|, 
cr = 0.40 dex, is indicated by a red dotted line. Bottom: Distribution of RMS scatter for a set of 10,000 simulated binary populations that 
were constructed by randomly pairing primaries with secondaries. As before, we show the RMS scatter of our observed population with a 
red dotted line; only 14 of the 10,000 simulated populations have tr<0.40 dex, indicating that our binary pairs are more coeval than Taurus 
at ~3<t significance. 

al. 1987); binary fragmentation is likely to occur dur- 
ing this collapse, so our limit is consistent with the pre- 
dicted formation timescale. Recent simulations that ex- 
ploit smoothed-particle hydrodynamic and N-body codes 
also suggest that most binary systems form within <0.5 
Myr (Delgado-Donate et al. 2004). 

In contrast, the timescale for global star formation is 
likely to be much longer, representing either the turbu- 
lent dissipation timescale (~1 Myr; Ballesteros-Paredes 
et al. 1999; Elmegreen 2000) or the ambipolar diffusion 
timescale (-3-10 Myr; Mouschovias 1976; Shu 1977). 
Our limit on the age dispersion of binary pairs is shorter 
than either timescale, while the overall age dispersion for 
unrelated pairs of stars is consistent with the ambipolar 
diffusion timescale, but only marginally with the turbu- 
lent dissipation timescale. We found a dispersion of ~0.6 



dex in |Alogr| for random pairs, corresponding to a dis- 
persion in logr of ~0.4 dex. For the median Taurus age 
of ~1.8 Myr, this corresponds to a typical age range of 
1-5 Myr. 

If our results do reveal two distributions (one popula- 
tion that appears coeval and one that does not), then 
the number of systems in each distribution will allow a 
constraint on the fraction of binary systems that are gen- 
uinely coeval. Of the 36 pairs of stars that we considered, 
24 are coeval to within |Alogr|<0.4 dex, while the other 
12 have ages which are more discrepant. This suggests 
that >67l 9 % of all binary systems are coeval with a dis- 
persion of <0.16 dex. However, many of the non-coeval 
pairs could be affected by systematic errors while being 
genuinely coeval, so this fraction is a lower limit. More 
intensive study of the apparently non-coeval pairs should 
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Fig. 4. — Difference in binary component age, |Alogr|, as a 
function of binary mass ratio. We see no evidence of a trend with 
q, as the standard deviation in |Alogr for q < 0.3 and for the full 
sample are the same for all pairs (0.41 dex versus 0.40 dex) and for 
the apparently coeval subset (0.17 versus 0.16 dex). 
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Fig. 5. — Difference in binary component age, |Alogr|, as a 
function of system mass. We see no trend for low-mass systems 
to appear more discrepant, but are unable to test systems with 
M >1.5 Mq and can only test a handful of systems with M <0.7 
M Q . 
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Fig. 6. — Difference in binary component age, |Alogr|, as a 
function of binary separation. The standard deviations in |Alogr| 
for the inner and outer halves are 0.37 dex versus 0.42 dex (for 
the full set) and 0.14 dex versus 0.18 dex (for the coeval subset); 
in both cases, the inner and outer halves are divided at 800 AU. 
This indicates that binary systems of all separations are similarly 
coeval. 

be a priority; as we discussed above, some stars (such as 
V710 Tau C, DK Tau A, and XZ Tau B) already seem 
potentially affected by systematic errors and might be 
rejected from our sample based on additional followup 
observations. 

Finally, our results indicate that the properties of a bi- 
nary system correlate only modestly with the formation 
timescale; even extreme systems (with very wide sepa- 



rations or disparate masses) appear similarly coeval on 
average. Conversely, unbound pairs of stars that are only 
modestly separated (~1 pc) show the full age dispersion 
of the association. These trends strongly indicate that 
binary coevality is a natural result of the binary forma- 
tion process itself, not a reflection of any trend for star 
formation to occur simultaneously within larger regions 
of the natal molecular cloud. 

6. SUMMARY 

We have studied the binary population of the Taurus- 
Auriga association in order to quantify the frequency and 
degree of noncoevality in young binary systems. After 
identifying and rejecting the systems that are known to 
be affected by systematic errors (such as further multi- 
plicity or obscuration by circumstellar material) , we used 
pre-main sequence evolutionary tracks to infer individual 
ages for the individual stars in binary systems, and hence 
the relative binary ages. We have found that the over- 
all dispersion in the relative ages (|Alogr|) is 0.40 dex, 
though the distribution actually appears bimodal. Ran- 
dom pairs of Taurus members are coeval only to within 
0.58 dex, suggesting that Taurus binaries are more coeval 
than the association as a whole. 

The bimodality indicates that our sample is comprised 
of two populations, with ^2/3 appearing to be coeval 
binaries with a dispersion of criAiogrl^O.16 dex and the 
other ~l/3 appearing to be systematically offset from 
coevality by ~0.6 dex. The non-coeval population shows 
no trends with respect to the system mass, separation, or 
mass ratio, which defies the predictions of many forma- 
tion scenarios for non-coeval systems. We therefore sug- 
gest that the non-coeval population is comprised mainly 
of unrecognized hierarchical multiples, stars seen in scat- 
tered light, or stars with NIR disk excesses; identifying 
any truly non-coeval systems will require additional fol- 
lowup to rule out or correct for these explanations. The 
full range of apparent ages in our sample is ~1 dex, which 
suggests that a binary system tends to form in a very 
short period of time relative to the global star formation 
timescale for Taurus. 

Finally, we found that the relative coevality of bi- 
nary systems does not depend significantly on the sys- 
tem mass, mass ratio, or separation. However, any pair 
of Taurus members wider than ~10' (~0.7 pc) shows the 
full age spread of the association. This suggests that the 
enhanced coevality is seen only for binary systems and 
not for neighboring stars that formed from separate pro- 
tostellar cores. The apparent coevality of a large fraction 
of our sample is also a partial endorsement of pre-main 
sequence isochrones. We did invoke several corrections 
to the model tracks, but any additional mass-dependent 
error in ages would cause systems with unequal mass ra- 
tios to appear less coeval, and we see no such trend to 
within the uncertainties in our results. 
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Fig. 7. — Difference in age as a function of (large-scale) separation for all possible pairs of Taurus members among our binary sample. We 
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of Taurus members (0.58 dex; Section 5.1). This result indicates that the coevality seen for binary systems is limited to smaller spatial 
scales, and perhaps only to binary systems themselves. 
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APPENDIX 
A. THE SINGLE STARS IN TAURUS 

The single stars of Taurus provide a useful check on the validity of our results, as well as providing their own 
constraints on its star formation history. In this appendix, we compile a sample of all stars which have a significant 
probability of being single (based on nondetections with one or more high-resolution imaging techniques). We then 
place these single stars on an HR diagram and estimate their ages and masses. Finally, we investigate the dependence 
of apparent age on stellar mass and on location within the association. 

A.l. Sample 

We list our sample of apparently single stars in Table 3, including all of the references for our adopted parameters 
(singleness, spectral type, and extinction). We based our sample on the compilation of all Taurus members that we 
originally described in Kraus & Hillenbrand 2007a), and then we omitted all stars that did not have at least one 
observation at high angular resolution. We then searched the literature for spectral types, requiring uncertainties of 
<1 subclass for spectral types >K0 and <2 subclasses for earlier-type stars. 

Most stars have only been surveyed for multiplicity to a separation limit of 50-100 mas (7-15 AU), so close binary 
systems still contaminate this sample. However, we do not expect any mass-dependent systematic biases. The binary 
fraction drops significantly with declining primary mass (e.g. Kraus et al. 2006, 2008), but much of that drop is 
seen among the wider binary systems (e.g. Kraus & Hillenbrand 2009) that would have fallen outside our required 
sensitivity limit of <15-20 AU. 

We inferred the fundamental properties of these stars (luminosity and temperature, then age and mass) using the 
methods described in Section 3. We list these properties in Table 3. 

A. 2. The Ages of Single Taurus Members 

In Figure Al, we show the HR diagram for our sample of single Taurus members. The composition of our single-star 
sample is significantly different from our binary sample, featuring many high-mass (>1 Mq) and low-mass (<0.3 M Q ) 
members, but few members with intermediate masses. This difference in composition is driven largely by selection 
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biases since few binaries at either extremum of mass have spatially resolved spectra. High-mass binaries in Taurus 
tend to be hierarchical multiples with additional components, while low-mass binaries were difficult to observe with 
spectroscopy before the recent advent of laser guide star AO. As a result, we must be very cautious in comparing the 
bulk properties of both samples. 

We also note that many members fall below our designated lower edge of the Taurus sequence (the 5 Myr isochrone) 
and might have erroneously low luminosities. Most of these members have not been well-studied (e.g. ITG33a and 
104301+2608), so we can not reject them with certainty, but their presence as extreme outliers in our plots invites 
skepticism. Most of the highest-mass members (SpT<K3; M >1 M ) also fall systematically below the 1-2 Myr 
isochrone. Few of our binary sample members fall in this mass range, but those that do are either obviously erroneous 
(the HBC stars) or appear genuinely young. We have no satisfactory explanation for this discrepancy between the 
single stars and binary components since an error in our methods or in the underlying models should affect both 
populations equally, but the small number of high-mass stars in our binary sample suggests that the single-star sample 
might provide a more reliable indication of the true empirical isochrone. 

In Figure A2, we show the model-derived age as a function of model-derived mass for our sample of single stars. 
Most stars seem to track the median age of Taurus, but as we noted above, the highest-mass stars (>1 M ) appear 
systematically older. The brown dwarfs of Taurus (M <0.07 M ) have very uncertain ages, so it is difficult to determine 
when they formed in relation to the stars. This uncertainty is driven by the physics of brown dwarf contraction, as 
isochrones at ages of <10 Myr follow similar tracks in the HR diagram. If we only consider members with masses of 
~0.07-0.9 M Q , then the median age of Taurus is log(r)= 6.25 ± 0.05 yr (1.8 ± 0.2 Myr). 

Finally, in Figure A3, we show the spatial distribution of our single-star sample on the sky, with the position of each 
star color-coded according to its age, as well as the mean age for the eastern subgroup, southern subgroup, and the 
eastern and western halves of the central filaments. The eastern subgroup appears ^4er older than the other subgroups 
(3.2 Myr, versus 1.4-1.9 Myr); if this age difference is genuine, then it suggests that star formation occurred first in 
Auriga, then in the rest of the association. However, this apparent age difference could also be a three-dimensional 
projection effect; the difference of ~0.25 dex in mean age could be explained if the distance to Auriga stars were 
~15-20% larger than the mean distance to Taurus (~170 pc versus ~145 pc). Otherwise, there is no significant trend 
in the ages of Taurus members, suggesting that global star formation proceeded nearly simultaneously (to within <0.3 
Myr). 

B. THE COEVALITY OF TRIPLE AND QUADRUPLE SYSTEMS 

High-order multiple systems are a critical tool for constraining stellar evolutionary models. If these multiple systems 
form coevally, then they provide a simultaneous test of the models at three or more masses. This feature was exploited 
by White et al. (1999) to constrain models with the well-known quadruple system GG Tau and to infer the best set 
of models to use for low-mass stars (the Lyon models) as well as to establish the best temperature scale for young 
stars (Luhman et al. 2003). We now extend this analysis to a quadruple system (UZ Tau), a quadruple that is part 
of a sextuple system (V955 Tau + LkHa332/G2), three components each of two systems that are not yet completely 
characterized (FV Tau and V773 Tau), and three triple systems (FS Tau, V710 Tau, and HL Tau/XZ Tau), plus we 
replicate the analysis of White et al. (1999) for GG Tau to provide context. 

As we show in Figure Bl, all three of the quadruple systems appear to have consistent ages. The consistency of GG 
Tau is partly a result of its previous role in calibrating stellar models and temperature scales, but UZ Tau appears to 
be almost coeval and similarly consistent. Three components of V955 Tau + LkHa332/G2 also fall along the 1 Myr 
isochrone, but V955 Tau B has an inferred age of ~3 Myr. This is ^2<r away from a consistent age, but among 12 
components, we would expect ~0.6 outliers at >2o\ 

Dynamical masses are available in the literature for the UZ Tau Aab and GG Tau Aab pairs and can be compared to 
those inferred from the HR diagram. The consistency is mixed. Guilloteau et al. (1999) found from the circumbinary 
disk kinematics that the total system mass for GG Tau Aab is 1.28±0.07 M ; the total mass predicted by theoretical 
models (1.37 M ) agrees to within ~7%. By contrast, Prato et al. (2002) reported dynamical masses for UZ Tau Aa 
and Ab of 1.02±0.06 M and 0.29±0.03 M , while the masses predicted by theoretical models are 0.61 and 0.30 M . 
The secondary mass agrees very well, but the discrepancy in the primary mass is very puzzling because its position 
in the HR diagram is virtually identical to that of GG Tau Ab, which has excellent consistency between observations 
and theory. Prato et al. explored the possible sources of this discrepancy in much greater detail, so we simply note its 
existence as proof that HR diagram analysis plays a critical, but incomplete role in constraining stellar evolutionary 
models. A full study of evolutionary models must include their dynamical masses (e.g. Schaefer et al. 2008), not just 
their temperatures and luminosities. Truly precise tests will also require direct measurement of radii (Stassun et al. 
2008) rather than indirect estimates from the Stefan-Boltzmann law and the observed luminosity and temperature. 

None of the three-component tests in our sample provide the same consistency seen among the quadruple systems, 
though as we described in Section 3, the two largest discrepancies are likely to be systematic. The edge-on disk host 
Haro 6-5B sits very far below the Taurus sequence, unlike FS Tau AB, while V710 Tau C might also be seen in 
scattered light. HL Tau is also seen in scattered light in the optical, but our inferred age based on its J magnitude 
seems consistent with that of XZ Tau A, suggesting that the central star of HL Tau might dominate its luminosity in 
the NIR. As we described in Section 2, XZ Tau B was suggested to be a possible binary pair by Carrasco-Gonzalez 
et al. (2009), which would explain its apparent overluminosity. FS Tau A and FS Tau B also appear moderately 
discrepant, sitting 1.5a on either side of the 1 Myr isochrone. FV Tau Aa and FV Tau Ba have very consistent ages, 
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Fig. Al. — HR diagram for all members of our single star sample. The dashed lines denote isochrones at 1 Myr (red), 5 Myr (green), 
10 Myr (blue), and 50 Myr (black). Most Taurus members fall along the 1-2 Myr isochrone, but many fall significantly below that level, 
perhaps due to the presence of an edge-on disk, undiscovered binary companion, erroneous observations. The highest- mass stars (>1 Mq) 
also fall systematically below the 1-2 Myr isochrone, suggesting either that the models might not be calibrated correctly in this regime or 
that these stars formed earlier in Taurus. 

but FV Tau Bb sits somewhat lower in the HR diagram; the only spatially resolved spectrum for FV Tau Bb is very 
noisy, so the apparent underluminosity could actually indicate that it has a later spectral type (M5.0-M5.5 rather than 
M3.5). 
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Fig. A2. — Age as a function of mass for all members of our single star sample. We also show the model-derived median age of Taurus 
(1.8 Myr; red line) as determined from our sample, plus representative error bars at four different masses (blue). The mass-dependent age 
of our sample tracks the overall median age except at the high-mass end (>1 Mq), where stars appear older, and at the low-mass end 
(<0.07 Mq), where the uncertainties become very large. 
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Fig. A3. — Spatial distribution for all members of our single star sample in the mass range that is well-calibrated (0.07-0.9 Mq), 
color-coded by age (red circles: <1 Myr; green squares: 1-3 Myr; blue triangles: >3 Myr). We also show the mean age for the eastern 
subgroup, southern subgroup, and the eastern and western halves of the core regions. The eastern subgroup appears ~4<r older than the 
other subgroups (3.2 Myr, versus 1.4-1.9 Myr) and contains no stars with an apparent age of <1 Myr, but otherwise there is no apparent 
pattern in the ages of Taurus members. 
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Binary Sample: Observed and Inferred Properties 
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24. 


.5 


2360±1 


K=8.74 


Ml±l 


0.8 


3705±145 


5.35 


Y? 


4, 


13, 22 


V927 Tau A 


4 


31 


23 


.82 


+24 


10 


52. 


9 


267±6.8 


K=9.31 


M3±0.5 


1.4 


3415±75 


6.00 


N 


8, 


9, 23 


V927 Tau B 


4 


31 


23 


.82 


+24 


10 


52. 


9 


267±6.8 


K=9.79 


M3.5±0.5 


0.9 


3340±75 


6.56 


N 


8, 


9, 23 


HL Tau 


4 


31 


38 


A A 

. 44 


+18 


13 


57. 


.7 


23310±70 


J=10.62 


K5±l 


7.4 


4350±265 


4.25 


Y 


10 


15 


XZ Tau A 


4 


31 


40 


.0 7 


+18 


1 Q 

lo 


57. 


2 


300.6±1.3 


K=8.36 


M2±l 


1.4 


3560±145 


4.98 


Y 


8, 


9 


XZ Tau B 


4 


Ol 


An 
4U. 


X) ( 


1 1 Q 


"1 Q 

lo 


5 / . 


2 


300.6±1.3 


K=7.80 


M3.5±0.5 


1.4 


3340±75 


4.52 


Y 


8, 9 


HK Tau A 


4 


31 


50. 


.57 


+24 


24 


18. 


1 


2342±61 


K=8.64 


Ml±0.5 


2.3 


3705±75 


5.08 


Y 


1, 


9, 10, 23 


HK Tau B 


4 


31 


50 


.57 


+24 


2'1 


18. 


1 


2342±61 


K=11.96 


Ml±0.5 


2.3 


3705±75 


8.40 


Y 


1, 


9, 10 


V710 lau A 


4 


31 


57. 




+ 18 


21 


QO 
OO. 


1 


O OO A 1 O 

3224±3 


K=9.38 


Ml±l 


0.9 


3705±145 


5.98 


Y 


1, 


13, 21 


V710 Tau B 


4 


31 


57. 




+18 


21 


38. 


1 


3224±3 


K=9.44 


M3±l 


0.9 


3415±145 


6.19 


N 


1, 13, 21 


V710 Tau C 


4 


31 


57 


.79 


+18 


21 


38. 


1 


28000±70 


J=12.26 


M3±0.5 


0.9 


3415±75 


8.00 


Y 


1, 3, 16 


GG Tau Aa 


4 


32 


30. 


.35 


+17 


31 


40. 


6 


10100±7 


J=9.07 


K7±l 


0.7 


4060±250 


4.50 


Y 


13, 17, 21 


GG Tau Ab 


4 


32 


30. 


.35 


+17 


31 


40. 


6 


250.2±2.6 


J=9.95 


M0.5±0.5 


3.2 


3775±75 


4.87 


Y 


13, 17, 21 


Tan Ra 


/[ 


32 


30 


.35 


+17 


31 


40. 


u 


1 m nn-t-7 


j = H 28 


1V1 U UZ u . o 


u.u 




7 94 


Y 


13, 14, 17, 21 


GG Tau Bb 


4 


32 


30 


35 


+17 


31 


40. 


6 


1476.5±6.5 


J=12.96 


M7.5±0.5 


0.0 


2795±85 


9.10 


Y 


13, 14, 17, 21 


UZ Tau Aa 


4 


32 


43 


.04 


+25 


52 


31. 


1 


SB 


H=8.54 


Ml±l 


1.5 


3705±145 


4.77 


Y 


1, 18, 22 


UZ Tau Ab 


4 


32 


43 


.04 


+25 


52 


31. 


1 


SB 


H=9.36 


M4±l 


1.5 


3270±145 


5.73 


Y 


1, 18, 22 


UZ Tau Ba 


4 


32 


43 


.04 


+25 


52 


31. 


1 


3539.5±2.1 


H=8.46 


M2±0.5 


0.6 


3560±75 


4.91 


Y 


8, 


9, 21 


UZ Tau Bb 


4 


32 


43 


04 


+25 


52 


31. 


1 


367.8±1 


H=9.18 


M3±0.5 


1.8 


3415±75 


5.48 


Y 


8, 


9, 21 


GH Tau A 


4 


33 


06. 


.22 


+24 


09 


34. 





311.1±1.3 


K=8.66 


M2±0.5 


0.0 


3560±75 


5.44 


Y 


8, 


9, 21 


GH Tau B 


4 


33 


06 


.22 


+24 


09 


34. 





311.1±1.3 


K=8.45 


M2±0.5 


0.5 


3560±75 


5.17 


Y 


8, 9, 21 


IS Tau A 


4 


33 


36 


79 


+26 


09 


49. 


2 


222.8±2.4 


K=8.82 


M0±0.5 


3.3 


3850±90 


5.08 


Y 


8, 


9 



TABLE 1 

Binary Sample: Observed and Inferred Properties 



IS Tau B 


4 


33 


36 


.79 


+26 


09 


49. 


2 


222.8+2.4 


K= 


=10.72 


M3. 5+0.5 


3 


.6 


3340+75 


7. 


19 


N 


8, 9 


HN Tau A 


4 


33 


39 


.35 


+17 


51 


52. 


.4 


3142+1 


K= 


=8.51 


K5+1 





.5 


4350+265 


4 


.82 


Y? 


1, 13, 20 


HN Tau B 


4 


33 


39 


.35 


+17 


51 


52. 


4 


3142+1 


K= 


=10.81 


M4.5+1 





.5 


3200+180 


7. 


.72 


Y? 


4, 13, 20 


IT Tau A 


4 


33 


54. 


70 


+26 


13 


27. 


5 


2416+8 


K= 


=8.12 


K3+1 


4 


.1 


4750+155 


3 


.84 


Y? 


4, 13, 22 


IT Tau B 


4 


33 


54 


.70 


+26 


13 


27. 


5 


2416+8 


K= 


=9.54 


M4+1 


4 


.1 


3270+145 


5 


.98 


Y? 


4, 13, 22 


Haro 6-28 A 


4 


35 


56 


.84 


+22 


54 


36. 





647+12 


K= 


=10.12 


M2+0.5 


2 


.3 


3560+75 


6 


64 


Y? 


8, 9, 24 


Haro 6-28 B 


4 


35 


56. 


.84 


+22 


54 


36. 


.0 


647+12 


K= 


=10.48 


M3.5+0.5 


1 


.9 


3340+75 


7. 


14 


Y? 


8, 9, 24 


2M04414565 


4 


41 


45 


.65 


+23 


01 


58. 





12370+70 


J= 


=10.74 


M3+0.5 





.0 


3415+75 


6. 


.72 




3, 16 


2M04414489 


4 


41 


44 


.89 


+23 


01 


51. 


3 


12370+70 


J= 


=14.42 


M8. 25+0.25 





.0 


2630+78 


10. 


56 




3, 16, 19 


LkHa332-G2 A 


4 


42 


07. 


.33 


+25 


23 


03. 


2 


234.1+4.5 


K= 


=8.38 


MO.5+0.5 


2 


.0 


3775+75 


4. 


.82 


N 


8, 9, 22 


LkHa332-G2 B 


4 


42 


07. 


33 


+25 


23 


03. 


2 


234.1+4.5 


K= 


=9.16 


M2. 5+0.5 


3 


.3 


3485+75 


5 


.61 


N 


8, 9, 22 


V955 Tau A 


4 


42 


07 


.77 


+25 


23 


11. 


8 


330.9+1.2 


K= 


=8.18 


K7+0.5 


2 


.8 


4060+125 


4. 


.26 


Y 


8, 9, 21 


V955 Tau B 


4 


42 


07 


77 


+25 


23 


11. 


.8 


330.9+1.2 


K= 


=9.72 


M2.5+0.5 


2 


.3 


3485+75 


6 


.28 


Y 


8, 9, 21 


UY Aur A 


4 


51 


47. 


.38 


+30 


47 


13. 


5 


878+17 


K= 


=7.68 


MO+0.5 





.6 


3850+90 


4. 


.24 


Y 


8, 9 


UY Aur B 


4 


51 


47. 


.38 


+30 


47 


13. 


5 


878+17 


K= 


=8.44 


M2.5+0.5 


2 


.7 


3485+75 


4. 


96 


Y 


8, 9 


2M04554757 


4 


55 


47 


.57 


+30 


28 


07. 


.7 


6310+70 


.1= 


=11.05 


M4.75+0.25 





.0 


3165+50 


7. 


15 




3, 7 


2M04554801 


4 


55 


48. 


01 


+30 


28 


05. 





6310+70 


J = 


=13.19 


M5.6+0.25 





.0 


3045+35 


9. 


31 




3, 7 


RW Aur A 


5 


07 


49. 


.54 


+30 


24 


05. 


1 


1417.5+3.4 


K= 


=7.25 


K2+2 


1 


.6 


4900+330 


3 


.11 


Y 


9, 10, 21 


RW Aur B 


5 


07 


49. 


.54 


+30 


24 


05. 


1 


1417.5+3.4 


K= 


=8.82 


K6+1 


1 


.6 


4350+265 


5 


.01 


Y 


9, 10, 21 



Note. — References: 1) Kenyon & Hartmann (1995), 2) Leinert et al. (1993), 3) Kraus et al. (in prep), 4) Duchcnc et al. (1999), 5) White et al. (in prep), 6) Boden et al. (2007), 7) Luhman 
(2004), 8) Hartigan & Kenyon (2003), 9) White & Ghez (2001), 10) White & Hillenbrand (2004), 11) Krist et al. (1998), 12) Konopacky et al. (2007), 13) Correia et al. (2006), 14) White & Basri 
(2003), 15) Krist et al. (1995), 16) Kraus & Hillenbrand (2009), 17) White et al. (1999), 18) Prato et al. (2002), 19) Luhman (2006), 20) Furlan et al. (2006), 21) McCabc et al. (2006), 22) 
Hartmann et al. (2005), 23) Luhman et al. (2006), 24) Andrews & Williams (2005). 



TABLE 2 
Derived Ages and Masses 



Name M (M Q ) log(r) (yr) 

(Adopted) 



M (M ) log(r) (yr) 
(Lyon) 



M (M ) log(r) (yr) 
(DM97) 



HBC 352 










>7.7 










>7.7 








>7.7 


fill/' Q F- Q 

HBC 353 










>7.7 










>7.7 








>7.7 


HBC 355 










>7.7 










>7.7 








>7.7 


HBC 354 










>7.7 










>7.7 








>7.7 


HBC 356 





79 


+0.08 
-0.03 


7 


68+° 
DS -0 


65 
27 










>7.7 





70+0. 09 
,y -0.05 


7 


56+° 


52 
27 


HBC 357 





79 


+0.08 
— 0.03 


7 


68+° 


65 
27 










>7.7 





7O+0.09 
,y -0.05 


7 


56+° 


52 
27 


V773 Tau Aa 




1. 


c+0.4 
-0.3 


5 


61 +0 
u -0 


22 
21 




1. 


3+0-3 
="-0.2 


5 92+ ' 28 

J ' az — 0.27 




1 5+°- 4 
l-°-0.3 


5 


61+° 
u -0 


22 
21 


V773 Tau Ab 




1. 


q+0.5 
u -0.1 


5 


62+2 


40 
24 




1. 


3+0-3 
-0.3 


6 25+ ' 48 
°- zd -0.31 




1 1+°- 4 

1 - 1 -0.4 


5 


691S 


33 
35 


2M04141188 


0.070j 


0.008 
0.007 


6 




19 
19 


0.0701 


0.008 
0.007 


6 46+ ' 19 


m +0.013 
U - 1±J --0.014 


6 


83tS 


11 

12 


FO Tau A 


0.4 


+0.10 
-0.17 


5 


861° 


21 

28 





40 


+ 
-0 


10 

17 


5 86+ ' 21 
o.oo_ 28 





1Q+003 
ly -0.03 


5 


20tS 


54 
27 


FO Tau B 


0.4 


+0.10 
— 0.17 


5 


86i° 


21 

28 





10 


+ 
-0 


10 

17 


5 86+°' oil, 

— 0.28 





1Q+0.03 
ly -0.03 


5 


20tS 


54 
27 


DD Tau A 





43 


+0.13 

— 0.23 


5 


631° 


23 
32 





43 


+ 
-0 


13 

23 


5.63+ ,' 23 

— (l-.iJ 





1O+0.03 
lo -0.02 


4 


85tS 


27 
25 


DD Tau B 





41 


+0.11 
— 0.19 


5 


791° 


22 
29 





41 


+ 
-0 


11 

19 


5.79+ ,' 22 

— 0.2)) 





1Q+003 
ly -0.02 


5 


091S 


35 
26 


FQ Tau A 





37 


+0.07 
— 0.05 


6 




23 
20 





37 


+ 
-0 


07 
05 


— 0.20 





nc:+0.05 
zo -0.03 


6 


2 ,+0 


20 
16 


FQ Tau B 





32 


+0.05 

— 0.05 


6 


431° 


20 
19 





32 


+ 
-0 


05 
05 


6.43+°' 2 ° 

— 0.19 





99 +0.03 
-0.03 


6 


2ltS 


17 
15 


LkCa 7 A 





67 


+0.04 
— 0.04 


6 


28±°o 


21 

21 





96 


+ 
-0 


09 
12 


6.61+ ;" 

— 0.22 





49+0.07 
* z -0.06 


5 


9ltS 


23 
18 


LkCa 7 B 





35 


+0.06 
— 0.07 


6 


17±8 


17 

22 





35 


+ 
-0 


06 
07 


6.17+"-" 

— 0.22 





91 +0.04 
zl -0.03 


5 


95tS 


15 
47 


FS Tau A 





66 


+0.03 

— 0.03 


5 


671° 


20 
20 


1 


06 


+ 
-0 


10 
12 


6.06+ ,' 2 . 2 

— 0.2o 





O9+0.05 
c>z -0.04 


5 


ietS 


42 
39 


FS Tau B 





32 


+0.05 
— 0.05 


6 


401° 


20 
19 





32 


+ 
-0 


05 
05 


6.40+ / 20 , 

—0.19 





99+0.03 
zz -0.03 


6 


istS 


17 
15 


Haro 6-5B 










>7.7 










>7.7 








>7.7 


FV Tau A 





83 


+0.06 

— 0.06 


5 




19 
18 


1 


54 


+0 
-0 


15 
15 


6 .38±S:i 





56 +0.l0 


5 


61±g 


20 
17 


FV Tau/c A 





54 


+0.05 
-0.07 


6 


081° 


19 
18 





55 


+ 
-0 


07 
07 


e.oslS;! 8 , 





9^+0.04 
zo -0.04 


5 


72tS 


16 
59 


FV Tau/c B 





30 


+0.07 

-0.05 


6 


671° 


26 
20 





30 


+ 
-0 


07 
05 


6.67+S: 2 2 6 o 





9O+0.04 
zo -0.03 


6 


48tS 


22 
19 


DF Tau A 





61 


+0.19 
-0.01 


5 




22 
24 





86 


+ 
-0 


06 
10 


c o 7 +0.17 
° °l -0.16 





1Q+0.02 
ly -0.02 


4 


22tS 


25 
24 


DF Tau B 





65 


+0.05 
-0.08 


5 


741° 


15 
27 





65 


+0 
-0 


08 
08 







99+0.03 
z -0.02 


4 


83tS 


28 
24 


2M04284263 A 





16 


+0.04 
-0.04 


6 


50±°, 


16 
19 





16 


+ 
-0 


04 
04 


6-50±g;ii 





1 7 +0.02 
1 ' -0.02 


6 


52tS 


20 
17 


2M04284263 B 


0.098j 


0.025 
0.019 


6 


621° 


34 
23 


o-098+S;SS 


6 62+ ' 34 

D ' D -0.23 


U - i4S -0.024 


6 


9ltS 


19 
18 


UX Tau A 




1. 


O+0.3 

J -0.4 


6 




26 
30 




1. 


3+ 


3.3 
3.3 


6 79+0- 19 
u - ,a -0.42 




1 O+0.3 
1 -°-0.4 


6 


iotS 


26 
30 


UX Tau C 





16 


+0.04 
-0.04 


6 


sitS 


16 

19 





16 


+0 
-0 


04 
04 


6 51 +016 





17+O.O2 
1 1 -0.02 


6 


54tS 


20 
17 


FX Tau A 





62 


+0.05 
-0.03 


5 


9ltS 


20 
19 





82 


+ 
-0 


19 
15 


6 14+0- 29 
D - 1 ^-0.27 





o n +0.07 


5 


54tS 


23 
59 


FX Tau B 





28 


+0.16 
-0.32 


5 


90+? 


33 
04 





28 


+0 
-0 


16 
32 


5 90+?' 33 

1.04 





10+0. 06 
lo -0.04 


5 


4ltS 


56 
38 


DK Tau A 





71 


+0.09 
— 0.06 


5 


34tS 


32 
36 


1 


30 


+0 
-0 


16 
22 


5.81+0'? 2 

— 0.80 





07+O.I3 
01 -0.12 


4 


76tS 


(i5 
51 


DK Tau B 





61 


+0.06 
— 0.04 


6 


17+° 
i '-o 


22 
20 





78 


+ 
-0 


18 
17 


6.38+ ,' 2 ? 

— 0.27 





oc+0.08 
OJ -0.07 


5 


sitS 


21 
23 


V927 Tau A 





42 


+0.07 
— 0.06 


6 


21±S 


19 
18 





42 


+ 
-0 


07 
00 


6.21+0-: 9 

— 0.18 





24 +0.04 
z -0.04 


5 


93tS 


15 
31 


V927 Tau B 





33 


+0.05 
— 05 


6 


32tS 


20 
18 





33 


+ 
-0 


05 

or, 


6.32±°?° 





99+0.03 
zz -0.03 


6 


12tS 


15 
15 


HL Tau 





83 


+0.11 

-0.10 


5 


93tS 


24 

25 


1 


53 


+ 
-0 


21 
27 


D - O1 -0.35 





c-7+0.21 
°'-0.15 


5 


68+S 


33 
22 


XZ Tau A 





59 


+0.07 
-0.04 


5 


setS 


20 
23 





67 


+ 
-0 


14 
12 


5 97+°' 27 
°- M '-0.28 





9C+O.O6 
zo -0.05 


5 


16tS 


57 
36 


XZ Tau B 





46 


+0.17 
-0.30 


5 


36tS 


24 
36 





46 


+0 
-0 


17 
30 







1 7+O.O2 
1( -0.01 


4 


44tS 


25 
25 


HK Tau A 





62 


+0.03 
-0.02 


6 


oil" 


19 
17 





80 


+ 
-0 


10 
09 


6 24+ ' 2 
D - z ^-0.21 





09+O.O5 
,3z -0.04 


5 


69+S 


15 
45 


HK Tau B 










>7.7 










>7.7 








>7.7 


V710 Tau A 





62 


+0.09 
-0.07 


6 


57+° 


24 
21 





73 


+0 
-0 


15 
18 


fi 74+ - 26 
D -'+-0.32 





40+O.12 
4J -0.10 


6 


2ltS 


31 

25 


V710 Tau B 





40 


+ 0.14 
-0.12 


6 


29IS 


32 
30 





40 


+ 
-0 


14 
12 


fi oq+0.32 
D - zy -0.30 





94+0.10 
z ^-0.06 


6 


oitS 


22 

35 


V710 Tau C 





35 


+0.07 
-0.07 


7 


24tS 


21 

25 





35 


+0 
-0 


07 
07 


- 24 +o.2i 

' - z4 -0.25 





9O+0.05 
zy -0.04 


7 


06tS 


24 
24 


GG Tau Aa 





73 


+0.09 
-0.08 


5 


88+S 


25 
30 


1 


21 


+0 
-0 


25 
22 


fi Q4+0.35 





,0 + 0.15 
+ J -0.12 


5 


621° 


25 
57 



TABLE 2 
Derived Ages and Masses 



GG Tau Ab 


64+ ' 03 
u - D4 -0.03 


5. 


93 


+0. 
-0 


19 
18 


0. 


91; 


fO.10 
-0.10 


6.23 + ° 


20 
22 


0. 


OO+0.05 
°°-0.04 


5 


.591° 


18 

46 


GG Tau Ba 


104+ ' 024 


5 


.80 


+0. 
-0 


31 
51 


o io4+°' 024 

u - lu4 -0.020 


5.80+° 


31 
51 


n ,,4+0.014 
U.lJ4_o.oi2 


6 


■nig: 


16 
44 


GG Tau Bb 


r)44 +0 ' 011 
u.U44_ () 013 


5. 


63 


+0. 
-2 


61 
.08 


044+ ' 011 
u - u + 4 -0.013 


5.631° 


.61 
08 


062+ 021 
U - UDZ -0.018 


6 


561° 


14 
71 


UZ Tau Aa 


n 62+ 05 


5 


.83 


+ 
-0 


19 

20 


0. 


.84; 


f0.18 
-0.14 


6.06+° 


2!) 
.27 


0. 


9O+0.07 
zy — 05 


5. 


.271° 


43 
42 


UZ Tau Ab 


o 26+ ' 20 


5. 


76 


+0. 

-1 


38 
20 


0. 


9 fi +0.20 
ZD -0.35 


5.76+? 


38 
20 


0. 


i 7 +0.06 
1 ' — 04 


5 


.201° 


67 
34 


UZ Tau Ba 


n 50+0 05 
u - oy -o.oi 


5. 


.81 


+0. 
-0 


19 
21 


0. 


f.n + 0.07 

b8 -0.08 


5.93+° 


19 
.17 





9^+0.04 
■^° — 03 


5 


091° 


53 
28 


UZ Tau Bb 


u -+ y -0.09 


5. 


98 


+0. 
-0 


18 
20 


0. 


49; 


f0.07 
-0.09 


5.98t° 


18 
.20 


0. 


99+0.04 
zz — 03 


5 


331° 


53 
31 


GH Tau A 


u - o5 -0.05 


6 


.15 


+0. 
-0 


19 
20 


0. 


6i; 


f 0.09 
-0.07 


6.20l° 


20 
.18 





9O+0.05 
zy — 04 


5 


.771° 


15 
47 


GH Tau B 


U - tl °-0.02 


5. 


98 


+0. 
-0 


20 
20 


0. 


64; 


f 0.08 
-0.07 


6.07+° 


18 
.18 


0. 


9 7 +0.04 
z '-0.03 


5 


.501° 


29 
49 


IS Tau A 


n 67+ 04 

u -°'-0.03 


6 


.11 


+0. 
-0 


21 
20 


0. 


99" 


f 0.09 
-0.12 


6.46l° 


19 
21 





oq+0.06 
■^-0.05 


5 


77+° 

— 


19 
18 


IS Tau B 


n 3n+ 006 


6. 


62 


+0. 
-0 


25 
20 


0. 


30; 


f 0.06 
-0.05 


6.62+°, 


25 
.20 


0. 


9O+0.04 

^•-"-o.os 


6 


.421° 


21 
18 


HN Tau A 


°- 85 -o'.io 


6 


.27 


+0. 
-0 


27 
25 


1. 


35" 


f 0.13 
-0.16 


6.85+° 


30 
34 





ec+0.23 
DO -0.18 


6 


051° 


38 
32 


HN Tau B 


u - zu -o.io 


6 


.59 


+0. 
-0 


45 
49 





20" 


f 0.13 
-0.10 


6.59t°, 


45 

49 


0. 


iq+0.07 
iy -0.05 


6 


■541° 


33 
23 


IT Tau A 


i o+°- 3 


5 


.96 


+0. 
-0 


32 
21 




1 S+ ' 2 
1 -°-0.2 


6.76+° 


18 
39 




1 0+°' 3 


5. 


.961° 


31 

32 


IT Tau B 


0.28±°; 3 5 2 


5 


.90 


+0. 

-1 


33 
04 





28" 


f 0.15 
-0.32 


5.901? 


33 
.04 





1O+0.06 
is -0.04 


5 


.421° 


55 
38 


Haro 6-28 A 


u - oz -0.08 


6 


.80 


+0. 
-0 


21 
22 





52" 


f 0.07 
-0.08 


6.80l° 


21 
22 





ocr+0.07 
" 3o -0.06 


6 


■4518 


24 
23 


Haro 6-28 B 


o-^O-o.os 


6 


.60 


+0. 
-0 


24 
20 


0. 


30" 


f 0.06 
-0.05 


6.60+° 


.24 
.20 





9O+0.04 
■^-O.OS 


6 


.391° 


21 
17 


2M04414565 


n o 7 +0.08 
U - d ' -0.05 


6 


.54 


+0. 
-0 


24 
20 





37" 


f 0.08 
-0.05 


6.541° 


24 
20 





9^+0.05 
zo -0.03 


6 


.261° 


21 
16 


2M04414489 


n27+ ' 006 
u - uz( -0.009 


6 


.47 


+0. 
-0 


44 
94 


027+ ' 006 
u - uz ' -0.009 


6.471° 


.44 
.94 


022+ ' 007 
u - uzz -0.006 




618 


82 
47 


LkHa332-G2 A 


n 64+ 03 

u - D4 -0.02 


5 


.90 


+0. 
-0 


19 

18 





92" 


+0.10 
-0.10 


6.20+° 


20 
21 





OO+0.04 
■ oo -0.0i 


5. 


.5618 


19 

50 


LkHa332-G2 B 


u - oz -0.07 


6 


.15 


+0. 
-0 


19 
18 


0. 


52" 


f 0.07 
-0.07 


6.151° 


19 
.19 





9f; +0.05 
ZD -0.04 


5 


•8118 


14 
53 


V955 Tau A 


n 79 +0.05 
u - '^-0.04 


5 


.73 


+0. 
-0 


22 
24 


1. 


24" 


+0.11 
-0.10 


6.20+° 


22 
22 





41 +0.07 
' 4i -0.06 


5 


51I8 


15 

55 


V955 Tau B 


n 4^+0.08 
u -+°-0.07 


6 


.46 


+0. 
-0 


23 
20 





45" 


f 0.08 
-0.07 


6.46+° 


23 
.20 





90+0. 06 
z °-0.04 


6 


•1218 


20 
16 


UY Aur A 


0.66I 002 


5 


.56 


+0. 
-0 


20 
20 


1. 


07" 


+0.12 
-0.12 


5.95+° 


23 
.25 





,1 +0.04 
• ol -0.04 


4. 


.941° 


50 
31 


UY Aur B 


u - DZ -0.08 


5 


.84 


+0. 
-0 


15 

23 





62" 


+0.08 
-0.07 


5.84+° 


19 
.20 





99+0.03 
zz -0.02 


4. 


.9918 


43 
24 


2M04554757 


n 20+ 03 

u - zu -0.03 


6 


.27 


+0. 
-0 


17 
17 


0. 


20" 


+0.03 
-0.03 


6.27l° 


17 
17 





1 7 +0.02 
0.01 


6 


.221° 


14 
13 


2M04554801 


n n92+ 008 


6 


.96 


+0. 
-0 


76 

27 


0.092I° ; o;« 


6.961°, 


76 
.27 


n ,00+0.015 

u.roo_o. 016 


7. 


1318 


16 
15 


RW Aur A 


i-4±8:? 


5 


.85 


+0. 
-0 


44 

.53 




n 1+0.3 
Z - L -0.5 


6.321°, 


54 
.42 




1 4+0-6 
i - 4 -0.7 


5. 


.8518 


44 
53 


RW Aur B 


0.86t° ;iJ 


6 


.40 


+0. 
-0 


26 
26 


1. 


2fi +0 .14 
ZD -0.14 


6.961°, 


29 
.33 





fiQ +0.24 
Dy -0.19 


6 


.1918 


39 

35 



TABLE 3 

Single Star Sample: Observed Properties 



Name 




RA 


Dec 

(J2000) 




J 

(mag) 


SpT 


Ay 
(mag) 


T eff 

(K) 


M bol 
(mag) 


M 
(-Me) 


log(r) 

(yr) 


Refs 


2MASSJ04080782+2807280 


4 


08 


07. 


82 


+28 


07 


28. 





12.44 


M3. 75+0.25 


1.0 


3305+35 


8.18 


0-25t°°I 


7.04«- 1 1 6 7 


1, 2 


LkCal 


4 


13 


14. 


14 


+28 


19 


10 


,8 


9.64 


M4+0.5 


0.0 


3270+70 


5.66 




5.70t° : g 


3, 4 


Anonl 


4 


13 


27. 


23 


+28 


16 


24. 


.8 


8.83 


M0+1 


3.6 


3850+178 


3.59 


0.65t° : {g 


5.13t° :2 2 ? 


3, 5, 6 


2MASSJ04141188+2811535 


4 


14 


11. 


88 


+28 


11 


53 


.5 


13.16 


M6.25+0.25 


0.7 


2960+28 


9.08 


0.069±g-88| 


6.50±g : ii 


1, 7 


FMTau 


4 


14 


13. 


.58 


+28 


12 


49. 


.2 


10.33 


M0+1 


1.9 


3850+178 


5.56 


0.68t°°t 


6.42t° : g 


3, 4, 6 


FNTau 


4 


14 


14. 


59 


+28 


27 


58. 


1 


9.47 


M5+0.5 


1.4 


3125+70 


5.23 


-0.20+° :2 ° 


4.201°^ 


3, 8 


CWTau 


4 


14 


17. 


00 


+28 


10 


57 


.8 


9.56 


K3+1 


1.9 


4730+155 


4.60 


i-oo±g-'g? 


6.39t° :2 ? 


3, 6, 8 


CIDA-1 


4 


14 


17. 


61 


+28 


06 


09 


,7 


11.73 


M5. 5+0.5 


3.0 


3060+68 


7.02 




4-40t 3 : «° 


9, 10, 11 


MHO-1 


4 


14 


26. 


40 


+28 


05 


59 


.7 


11.52 


M2.5+0.5 


5.7 


3490+73 


5.90 


0.49t° : ^ 


6.30±° : ?° 


6, 10, 12 


FPTau 


4 


14 


47. 


31 


+26 


46 


26. 


,4 


9.89 


M4+0.5 


0.2 


3270+70 


5.85 


0.27t° : ^ 


t- oo+0.26 
0-8^_0 31 


3, 4 


CXTau 


4 


14 


47. 


86 


+26 


48 


11 


.0 


9.87 


M2.5+0.5 


0.8 


3490+73 


5.59 


n p-o+0.05 
u *°° — 07 


6-15t° : i! 


3, 4 


KPNO-Tau-1 


4 


15 


14. 


71 


+28 


00 


09 


,6 


15.10 


M8.5+0.25 


0.4 


2555+78 


11.13 


023+ 010 

U.UiiO _ Q QQQ 


6.75t° : *° 


13, 14, 15 


2MASSJ04152409+2910434 


4 


15 


24. 


09 


+29 


10 


43 


.4 


13.68 


M7+0.25 


2.0 


2880+33 


9.27 


n Qro+0.006 
u.uuo_q 005 


6 S1+ 26 
D.01_o 26 


1, 2 


2MASSJ04161210+2756386 


4 


16 


12. 


10 


+27 


56 


38 


.6 


12.27 


M4.75+0.25 


2.0 


3160+38 


7.81 


n 1S1 +0.022 

u.iei_ 


6.55_o 14 


7, 16 


2MASSJ04161885+2752155 


4 


16 


18. 


85 


+27 


52 


15. 


.5 


12.55 


M6.25+0.25 


1.0 


2960+28 


8.39 


0.080±° : 8g| 


e.istg'lS 


1, 2 


2MASSJ04163911+2858491 


4 


16 


39. 


11 


+28 


58 


49. 


1 


12.72 


M5. 5+0.25 


2.8 


3060+38 


8.07 


0-13tS-g§ 


6-44t° :i 1 i 


1, 2 


CYTau 


4 


17 


33. 


73 


+28 


20 


46 


.9 


9.83 


Ml+0.5 


0.1 


3705+73 


5.64 


n C1 +0.04 

0.61I 02 


6.35t° :2 ° 


3, 8 


KPNO-Tau-10 


4 


17 


49. 


55 


+28 


13 


31 


.9 


11.89 


M5+0.25 


0.0 


3125+33 


8.02 


15+ 002 


6.55t° :i 13 


1, 14, 17 


V410-Xrayl 


4 


17 


49. 


65 


+28 


29 


36 


3 


11.02 


M3. 75+0.25 


0.9 


3305+35 


6.78 


0.29ioT 2 




6, 10, 18 


V410-Anonl3 


4 


18 


17. 


11 


+28 


28 


41. 


9 


12.96 


M6+0.5 


3.8 


2990+63 


8.02 


0.09t° : g 3 2 




6, 12, 15 


KPNO-Tau-11 


4 


18 


30. 


31 


+27 


43 


20 


.8 


11.89 


M5.5+0.25 


0.0 


3060+38 


8.01 


0.13±° : °I 


6.42t°' 1 i 2 7 


1, 17 


KPNO-Tau-2 


4 


18 


51. 


16 


+28 


14 


33 


,2 


13.92 


M7.5+0.25 


0.4 


2795+45 


9.96 


0.040±g-ggl 


6 71 + 013 
u - ' x — 67 


13, 14, 15 


HBC376 


4 


18 


51. 


70 


+17 


23 


16 


.6 


10.03 


K7+1 


0.0 


4060+250 


5.66 


0.79t°oi° 


a CQ +0.'27 
6 - 63 -0 26 


3, 4 


104158+2805 


4 


18 


58. 


14 


+28 


12 


23 


,5 


13.78 


M6+1 


8.6 


2990+123 


7.50 


in+ 07 

U.1U_Q 04 


5 91+ 45 


10, 19 


KPNO-Tau-12 


4 


19 


01. 


27 


+28 


02 


48 


,7 


16.31 


M9+0.25 


0.5 


2400+75 


12.31 


0-032±g-g?S 


7.46t° :2 l 


15, 17 


V410-Xray5a 


4 


19 


01. 


98 


+28 


22 


33 


,2 


11.99 


M5.5+0.5 


2.6 


3060+68 


7.41 




6.20t° : ^ 


6, 10, 12 


BPTau 


4 


19 


15. 


84 


+29 


06 


26 


.9 


9.10 


K7+1 


0.5 


4060+250 


4.59 


r, 7 o+0.09 
U- i <J_o 08 


5.94t° :2 | 


3, 5 


2MASSJ04202555+2700355 


4 


20 


25. 


.55 


+27 


00 


35 


.5 


12.86 


M5.25+0.25 


2.0 


3095+33 


8.42 


n 19 o+0.0l8 
— 014 


6.62t° : ii 


1, 7 


J2-157 


4 


20 


52. 


73 


+17 


46 


41 


.5 


11.62 


M5.5+0.5 


0.0 


3060+68 


7.74 


0.14±g-°i 


6-32+°^ 


3, 10 


CFHT-Tau-19 


4 


21 


07. 


95 


+27 


02 


20 


4 


13.85 


M5.25+0.25 


7.3 


3095+33 


7.95 


0.145t° ° 2 ° 


6.4610-11 


14, 16 


2MASSJ04213460+2701388 


4 


21 


34. 


60 


+27 


01 


38. 


.8 


11.90 


M5.5+0.25 


1.8 


3060+38 


7.53 


o.ut°-° t 


6.24t° : » 


7, 16 


CFHT-Tau-10 


4 


21 


46. 


.31 


+26 


59 


29 


6 


13.82 


M5. 75+0.25 


2.0 


3020+35 


9.38 


0.084t°-° 1 1 2 2 


6.88t°'g 


1, 2 


2MASSJ04215450+2652315 


4 


21 


54. 


50 


+26 


52 


31. 


5 


15.54 


M8.5+0.25 


1.0 


2555+78 


11.40 


0.027t° : °^ 


7-09t° : f 5 


1, 2 


DETau 


4 


21 


55. 


64 


+27 


55 


06. 


1 


9.18 


M2+0.5 


0.6 


3560+70 


4.93 


0.59±g-g* 


c- S9 +0.19 
— 20 


3, 8 


RYTau 


4 


21 


57. 


40 


+28 


26 


35 


.5 


7.16 


Kl+1 


1.8 


5080+175 


2.09 


o 1+0.6 
z - l -0.6 


5 61 + '' 28 


3, 8 


HD283572 


4 


21 


58. 


84 


+28 


18 


06. 


,6 


7.42 


G5+2 


0.4 


5770+98 


2.54 


2 0+ ' 4 
^• u -0.2 


6.62t°; 2 ° 5 


3, 5 


CFHT-Tau-14 


4 


22 


16. 


44 


+25 


49 


11 


.8 


13.06 


M7.75+0.25 


0.5 


2750+43 


9.06 


n nQC+0.005 
°- 036 -0.OO6 


4.l9t 2 ;l 2 


1, 2, 14 


CFHT-Tau-21 


4 


22 


16. 


76 


+26 


54 


57. 


1 


11.58 


Ml.5+0.25 


3.0 


3630+38 


6.63 


^8+ 02 
u - o5 -0.03 


6 Q1+ 018 
D ' yi -0.18 


2, 16 


2MJ04230607+2801194 


4 


23 


06. 


07 


+28 


01 


19. 


.4 


12.24 


M6.25+0.25 


0.0 


2960+28 


8.36 


n nsn+ 008 
O.OSO_ oos 


6 17+ 016 
D - 1 '-0.16 


1, 20 


CFHT-Tau-9 


4 


24 


26. 


46 


+26 


49 


50 


4 


12.88 


M5. 75+0.25 


0.5 


3020+35 


8.86 


0.087t°;^ 


6 57+ 21 
°-°'-0.17 


1, 2 


IPTau 


4 


24 


57. 


08 


+27 


11 


56. 


.5 


9.78 


M0+1 


0.2 


3850+178 


5.47 


0-68t°;°l 


c 0^7+0.23 
u -°' -0.24 


3, 5 


KPNO-Tau-3 


4 


26 


29. 


39 


+26 


24 


13. 


.8 


13.32 


M6+0.25 


1.6 


2990+30 


8.98 


0.077i°;°°t 


6 55+° 17 


13, 15 



Single Star 



TABLE 3 
Sample: Obsi 



KPNO-Tau-13 


4 


26 


57. 


33 


+26 


06 


28. 


.4 


11.28 


M5+0.25 


HBC388 


4 


27 


10. 


56 


+17 


50 


42 


.6 


8.79 


Kl±l 


KPNO-Tau-4 


4 


27 


28. 


00 


+26 


12 


05. 


3 


15.00 


M9. 5+0.25 


2MASSJ04290068+2755033 


4 


29 


00 


68 


+27 


55 


03. 


3 


14.02 


M8.25+0.25 


KPNO-Tau-5 


4 


29 


45. 


68 


+26 


30 


46. 


.8 


12.64 


M7.5+0.25 


IQTau 


4 


29 


51. 


56 


+26 


06 


44 


.9 


9.42 


MO.5+0.5 


CFHT-Tau-20 


4 


29 


59. 


.51 


+24 


33 


07 


9 


11.68 


M5+0.25 


KPNO-Tau-6 


4 


30 


07. 


24 


+26 


08 


20 


.8 


15.00 


M8.5+0.25 


CFHT-Tau-16 


4 


30 


23. 


65 


+23 


59 


13. 





14.96 


M8.25+0.25 


KPNO-Tau-7 


4 


30 


57. 


19 


+25 


56 


39 


.5 


14.52 


M8.25+0.25 


JH56 


4 


31 


14. 


44 


+27 


10 


18. 





9.70 


MO.5+0.5 


MHO-9 


4 


31 


15. 


.78 


+18 


20 


07 


.2 


11.21 


M5+0.5 


2MASSJ04311907+2335047 


4 


31 


19. 


07 


+23 


35 


04 


.7 


13.51 


M7.75+0.25 


MHO-4 


4 


31 


24. 


.06 


+18 


00 


21 


.5 


11.66 


M7+0.5 


CFHT-Tau-13 


4 


31 


26. 


69 


+27 


03 


18 


.8 


14.83 


M7.5+0.25 


LkHa358 


4 


31 


36. 


13 


+18 


13 


43 


3 


12.80 


M5.5+0.5 


HLTau 


4 


31 


38. 


44 


+18 


13 


57. 


,7 


10.62 


K5+1 


Jl-665 


4 


31 


58. 


44 


+25 


43 


29 


.9 


10.59 


M5+0.5 


2MASSJ04320329+2528078 


4 


32 


03. 


29 


+25 


28 


07 


.8 


11.72 


M6.25+0.25 


L1551-51 


4 


32 


09. 


27 


+17 


57 


22 


.8 


9.70 


K7+1 


Haro6-13 


4 


32 


15. 


41 


+24 


28 


59. 


,7 


11.24 


MO+0.5 


MHO-5 


4 


32 


16. 


.07 


+18 


12 


46 


.4 


11.07 


M7+0.5 


CFHT-Tau-7 


4 


32 


17. 


86 


+24 


22 


15. 





11.54 


M5. 75+0.25 


MHO-6 


4 


32 


22. 


11 


+18 


27 


42. 


.6 


11.71 


M4.75+0.25 


2MASSJ04322329+2403013 


4 


32 


23. 


29 


+24 


03 


01 


3 


12.34 


M7.75+0.25 


MHO-7 


4 


32 


26. 


.28 


+18 


27 


52 


.1 


11.11 


M5.25+0.25 
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NOTE. — References: 1) Kraus et al. (in prep), 2) Luhman (2006), 3) Kenyon & Hartmann (1995), 4) Leinert et al. (1993), 5) Tanner ct al. (2007), 6) Luhman et al. (2000), 7) Luhman (2004), 
8) Ghez et al. (1993), 9) White & Basri (2003), 10) White ct al. (in prep), 11) Hcrczcg fc Hillenbrand (2008), 12) Bricciio ct al. (1998), 13) Briccno et al. (2002), 14) Guicu ct al. (2006), 15) 
Kraus ct al. (2006), 16) Konopacky et al. (2007), 17) Luhman ct al. (2003), 18) Strom & Strom (1994), 19) White & Hillenbrand (2004), 20) Luhman et al. (2006), 21) Hartigan et al. (1994), 22) 
Sartoretti et al. (1998), 23) Hartigan ct al. (1995), 24) Simon et al. (1995), 25) Martin et al. (2001), 26) Slcsnick ct al. (2006), 27) Briccno ct al. (1999). 
a Some systems which sit extremely high or low in the HR diagram exceed the limits of the models and have nonphysical derived quantities (i.e. negative masses). 
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